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ABSTRACT 

Recently, the Kepler spacecraft has detected a sizable aggregate of objects, characterized by giant- 
planet-like radii and modest levels of stellar irradiation. With the exception of a handful of objects, 
the physical nature, and specifically the average densities, of these bodies remain unknown. Here, we 
propose that the detected giant planet radii may partially belong to planets somewhat less massive 
than Uranus and Neptune. Accordingly, in this work, we seek to identify a physically sound upper 
limit to planetary radii at low masses and moderate equilibrium temperatures. As a guiding example, 
we analyze the interior structure of the Neptune-mass planet Kepler-Z^d and show that it is acutely 
deficient in heavy elements, especially compared with its solar system counterparts. Subsequently, we 
perform numerical simulations of planetary thermal evolution and in agreement with previous studies, 
show that generally, 10 — 2OM0, multi-billion year old planets, composed of high density cores and 
extended H/He envelopes can have radii that firmly reside in the giant planet range. We subject 
our results to stability criteria based on extreme ultraviolet radiation, as well as Roche-lobe overfiow 
driven mass- loss and construct mass-radius relationships for the considered objects. We conclude by 
discussing observational avenues that may be used to confirm or repudiate the existence of putative 
low mass, gas-dominated planets. 
Subject headings: planets and satellites: interiors, planets and satellites: physical evolution 



1. INTRODUCTION 

The ever-growing transit data set collected by the 
Kepler spacecraft has proven to be instrumental to the 
advancement of our understanding of the properties of 
planetary systems. Thanks to the sheer size of the data 
set (~ 2500 planetary candidates as of Quarter 6) and 
the associated statistical ability to determi ne the charac- 
teristics of typ ical planetary systems ( jHoward et al.|201Q 



|Youdin||2011|), as well as hi^hli^ht some unex ptected ex- 
amples (l^oyl e et"ar]|2011 Welsh et al.||2012 ), important 
insights into plane t formation have already bee n gleaned 
from the analysis (Wolfgang & Laughlin 2011). 

It is interesting to note that the Kepler data set con- 
tains objects whose radii are similar to that of Jupiter 
(and in a few cases even exceed it substantially), in 
the moderate irradiation range (200K < Teq < 800K) 
(see Figure IT]). Specifically, the latest application of 
the pipeline to the sample suggests that of 1333 total 
planetary candidates in this T^.^ range, 68 have radii 
in the (i^sAT < R< 2i?jup) 



eq range 

range and 25 have radii 



th at exc eed Rj by more than a facto r of 2 (Batalha et 
aLl|2012|). Although the analysis of pemory ^ Seager 



suggests that a dominant portion ol the exces- 
sively large objects in the Kepler inventory are false pos- 
itives, the physical nature of objects characterized by 
Jupiter-like radii is of considerable interest. Neverthe- 
less, even basic information such as the average density 
is difficult to acquire since the overwhelming majority 
of the stars in the Kepler field are rather faint, making 
radial-velocity follow up observationally expensive. Bar- 
ring (near-) resonant syst ems, where transit timing vari- 
ations can be significant ( Holman fc Murray||2005 ), this 
means that the masses of the planets within the Kepler 



sample will remain observationally unconstrained and 
theoretical inquiries are desirable. 

It is well known that giant planets comprising hun- 
dreds of Earth masses can have large radii, that exhibit 
only weak dependent on mass, and are instead primarily 
controlled by their chemical composition and the inte- 
rior thermal state (Zapolsky & Salpeter 1969; Stevenson 
1982a). Furthermore, it is firmly established that radii 
of gaseous planets can increase with decreasing mass, 
thanks to the associated softening of the equation of state 
(Stevenson 1982a). Although, as illustrated by the wide- 



ranging numerical calculations of [Fortney et al.| (|2007|) , 
whether the radius increases or decreases with niass and 
the extent to which it does so, are rather sensitive to the 
amount of irradiation received by the planet as well as 
its chemical composition. 

In extreme proximity to the host star, the upturn in ra- 
dius is well pronounced. For example, an evolved 2OM0 
planet irradiated at Teq !::^ 1300K is roughly twice as large 
as its isolated counterpart (Baraffe et al. 2008). Depend- 
ing on the planetary age, at even higher temperatures 
(e.g. Teq = 200 0K), the discre pancy may be as large as a 



otlp 



factor of a few (Guillot 2005). On the other hand, plan- 



etary radii at iea 5-, lOOK, do not differ from those of 



isolated objects much ( [Fortney et al.|2007 [). The current 
observational frontier lies in between these extremes, and 
to da te, with the ex c eption of only a handful of studies 



(e.g. Rogers et al. 
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( 2011[) ), this parameter regime re- 
mains largel3nme5cplore37ln this study, we shall perform 
calculations that will place meaningful constraints on the 
mass-radius relationship of sub-Saturnian objects in the 
moderate irradiation regime. Specifically, the identifica- 
tion of a physically sound upper limit to the planetary 
radii at low masses and moderate equilibrium tempera- 
tures is the primary aim of this study. 
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Fig. 1. — Planetary radii as a function of planetary equilibrium ir- 
radiation temperature in the Kepler sample. Note the considerable 
presence of giant-planet-like radii in this irradiation regime. The 
data was obtained from http://planetquest.jpl.nasa.gov/kepler 

The possible range of chemical compositions of plan- 
ets is generally not well known. However, the relatively 
low densities exhibited by some members of the well- 
characterized subset of the Kepler catalog suggest that 
low overall metallicities cannot be ruled out. A par- 
ticularly important exampl e is the planet Kepler-SOd 
(Sanchis-Ojeda et al. 2012) which as we show below, 
has an envelope whose density does not exceed that of 
a cosmic H/He mixture substantially and cannot pos- 
sess a core as massive as that typically i nvoked in the 
core- accretion model of planet formation ( [Pollack et al. 
1996). Thus, motivated by the inferred structure ot Ke- 
p/er-30d, we shall limit ourselves to a consideration of 
the most favorable planetary compositions for the fabri- 
cation of large radii. That is, for definiteness and sim- 
plicity, here we focus on planets with well-defined cores 
and H/He gaseous envelopes, though it is possible that 
in real objects the core material is parti ally mixed into 
the envelop^ ( [Leconte fc Chabrier|[2QT2| . 

The paper is organized as tollows. In section 2, we de- 
scribe the setup of our numerical experiments and per- 
form simulations of planetary thermal evolution to ex- 
plore the interior structure of Kepler-SOd. In section 
3, we extend our calculations to lower masses and con- 
struct generic mass radius relationships, constrained by 
the hydrodynamical stability of the considered planets. 
We conclude and discuss our results in section 4. 

2. THE STRUCTURE OF KEPLER-30B 



Following initial detection (Batalha et al. 2012), the 
Kepler-30 system was st udied in greater detail by 
Sanchis-Ojeda et al. ('2012'), who determined the mass, 
radius, equilibrium temperature and age of Kepler-SOd 
to be M = 23.1 ± 2.7Me, R = 8.8± 0.5i?®, Teq = 364K 
(assuming zero albedo) and 2.0 ± 0.8 Gyr respectively. 
Here, we shall adopt the observed best fit parameters at 
face value for the generation of interior models. 

Naturally, any model we consider is subject to hydro- 
static equilibrium. With the knowledge of the equation 

^ Mixing of core and envelope may not change the radius much 
for Jupiter mass planets when they are compared at similar tem- 
peratures. However, the consequences of this mixing are in general 
not simple for the radius-mass relationship because it affects the 
cooling history of the planet as well as the density distribution for 
a given temperature. 
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Fig. 2. — The radius-core mass relationship for Kepler-30d, as- 
suming various envelope compositions. Blue curves correspond to 
envelopes with Z = 0.02 while red curves correspond to Z = 0.04. 
The cartoon in the left-bottom corner of the figure represents the 
considered two-layer interior models (here drawn to scale with a 
5M0 core). 

of state and an assumed radiative structure of the at- 
mosphere, the construction of a static interior model is 
relatively straight forward. This is however not enough, 
since the thermal state of the planet chan ges in time du e 
to radiative losses of the interior entropy ( Guillot|[l999 ). 
By extension, the planetary radius also contracts. TEus, 
in order to obtain definitive results that are characteris- 
tic of multi-Gyr old planets, evolutionary calculations of 
planetary structure are required. 

For our numerical experiments, we utilized the ME SA 
stellar and planet ary evolutio n software packag e ( P axton l 
et al||2Qll[[2013D . Following lBodenhei mereFal.| ( |200T| 7 
all ot our models comprised constant density [Pcore = 5 
g/cc) solid cores embedded in gaseous H/He envelopes. 
The baseline heat-flux arising from radioactive decay 
within the cores was taken to be 10~^ ergs/s/g, similar 
to that of the Earth. The envelope metallicity was varied 
between Z = 0.02 and Z = 0.04, while solar Y = 0.27 
and slightly super-solar Y = 0.35 values of the He mass 
fraction were explored. We note that if hydrodynamic 
mass-loss play ed a significant role in shaping the plane- 
tary structure ( Owen fc Wu|2013 ), a super-solar value of 
Y can in principle originate troni a preferential blow-off 
of Hydrogen. 

The analytical radiative equilibrium model o f the outer 
atmosphere was adopted from the work of (Guillot fc 
Havel||2011|). In the radia tive portion of the atmosphere 
followmg Guillot (2010) we take the ^ — --^-^-^ ^-^^ui^ 
and infrared opacities to be hzy = 10 " cm 
^IR = 4 X 10~^ cm^ g~^ respectively. These 
yield the closest agreement between the analytical radia- 
tive model used here and the state of the art numeri- 

( [2QQ8t . 
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choices 



cal models of 
much greater 



Fortney et al.| ([2008|. At optical depths 
han unity, the tabulated Rosseland mean 



opacities of Freedman et al.| ([2008) were used and the 



radiative convective boundary was computed as dictated 
by the Schwartzchild criterion. The reported radius of a 
given planet was taken to be the value corresponding to 
a chord optical depth of unity in visible light. 

The search for admissible models of Kepler-SOd was 
performed in the following way. For a given choice of Y 
and Z, the core- mass was varied between Mcore = ^Mm 
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and Mcore = 23M0. The resulting initial conditions were 
integrated forward in time, yielding a sequence of model 
radii that decrease monotonically with Mcore- Impor- 
tantly, radii also decrease monotonically with enhanced 
mean molecular weight, which means that there exists 
a maximum value of Mcore above which the planetary 
radius cannot be matched. There also exists a maximal 
extent to which the mean molecular weight of the enve- 
lope can exceed that of a cosmic H/He mixture. However, 
such core-less solutions are strongly disfavored because 
the mass of Kepler-SOd is too low for formation by grav- 
itational instability to be plausible. 

The R — Mcore sequences for various compositions 
are shown in Figure ([2|. Adopting a solar composi- 
tion envelope yields an upper bound on the core-mass of 
Mcore — 7M0. Meanwhile, the corresponding value for a 
Y = 0.35, Z = 0.04 envelope is a mere Mcore — S.SMq. 
Unfortunately, because the planet's gravitational har- 
monics are not known, no useful lower bound on Mcore 
can be formulated. It is noteworthy however, that the 
upper bound on Mcore is surprisingly low. 

The dominantly gaseous interior structure we obtained 
for Kepler-SOd is in sharp contrast with the inferred 
heavy elemen t-dominated interior structures of Neptune 
and Uranus (Fortney et al. 2011). This suggests that 
the diversity m composition and overall interior struc- 
ture of low-mass planets is generally much more exten- 
sive than what is captured within modern state-of-the-art 
core accretion models. More specifically, this implies that 
the nucleated instability mechanism can operate even for 
comparatively small cores. 

3. GENERIC MASS-RADIUS RELATIONSHIPS 

Motivated by the results attained above, in this sec- 
tion we construct generic mass-radius relationships for 
evolved planets with a specified composition, extending 
down to minimum feasible masses. The radius-mass se- 
quences were generated in a similar manner to the numer- 
ical experiments reported in the previous section. How- 
ever, with the aim to constrain the planetary radii from 
above, the compositions of the envelopes were kept solar 
{X = 0.71,7" = 0.27, Z = 0.02) across the models. The 
core mass was varied between 1,3 and 5 Mq, while the 
total planetary mass range of up to 0.1 Mjup (~ 33M0) 
was explored. 

Conventional generation of initial conditions within the 
framework of thermal evolution calculations is known 
to encounter numerical instabilities at sufficiently low 
masses. Consequently, here the initial conditions were 
constructed by imposing a slow mass-loss on a M = 
O.lMjup model. After the desired mass was attained, 
we imposed energy dissipation to the core and re-heated 
the gaseous envelope to the point where the thermal 
and gravitational energies of the body are comparable. 
The duration of the evolutionary sequences was formally 
taken to be 5 Gyr. However, it should be noted that 
the changes in planetary structure were relatively small 
after the first ~Gyr of integration. Likewise, we found 
the evolved radii to be largely independent of the de- 
tailed state of the i nitial condition, in agreement with 
publis hed literature ( Bodenheimer et al.||2001 Hubickyj 
eTaLJIMosl ). 

Wot all generated planetary models are guaranteed to 
be long-term stable. Indeed, some of the models we con- 



structed were characterized by radii, exceeding that of 
Jupiter by as much as a factor of a few, rendering their 
stability against evaporation questionable. Accordingly, 
we formulated a criterion for model rejection in terms of 
the mass-loss rate due to atmospheric escape. 

Irradiated extrasolar planets can be susceptible to 
mass- loss due hydro dynamic wind s originatin g in the 



mg 
3E 
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upper atmosphere ( . Murray-Clay et al. 

et al. 2010 [Lopez et al. ' 2012[ ). Such winds are gener- 
ated throuffn the photoiomzation of H (and the associ- 



photoiomzation of H (and 
ated heating) by extreme ultraviolet radiation. Provided 
that downward conductive heat fiux or radiative coolin g 
by H3" is not overwhelming (Murray-Clay et al. 2009), 
the characteristic timescale for energy- Im iited evapora- 
tion is given to an orde r of magnitude by (Watson et al. 
TMj [Yelle eni][2QQ8l ) 



eTrFEuv^luv 



(1) 



where G is the gravitational constant, Feuv = 
4.1(a/lAU)~^ erg s~^ cm~^ is the typi cal extreme ultra- 



violet fiux of a 5 Gyr old Sun-like star ( Ribas et al.|2005" 
Sanz-Forcada et al.J 2010'l ), Ktide = l-3(ii:Euv/^Hiiij/2 



(i?EUv/^Hiii)^/2 is a geometrical factor that accounts for 
the fact any given parcel o f gas only needs to reach the 
Hill radius to escape (Erkaev et al. 2007), and Revv 
is a radius at which the atmosphere becomes optically 
thick to extreme ultraviolet radiation i.e. nHaxw ^ 1 
where n is the atmospheric number density, T-L is the 
scale- height, and axuv — 10~^^ c m^ is the photoion- 
ization cross-section for Hydrogen (Murray-Clay et al. 
2009[ ). Meanwhile, e o:^ 0.25 is a factor that parameter- 
izes the efficiency of atmospheric escape. 

The EUV fiux is considered to be constant here since 
we are not seeking to model loss and stability during 
early epochs of evolution. However, any model that 
we deem stable at t = 5 Gyr will also likely be sta- 
ble at any time greatly exceeding the T-Tauri phase of 
the evolutionary sequence (e.g. t > 100 Myr), because 
we generally find characteristic loss timescales of order 
'?"e-iim ^ 100 Gyr or greater. 

In all our models, Rew never exceeded the exobase 
(a radius at which the molecular mean free path be- 
comes comparable to H), meaning that the atmospheres 
were never truncated by Jean's escape. However, for cer- 
tain models, Rew exceeded i^Hiii, implying mass-loss 
by Roche- lobe overfiow. In such ca ses the characteristic 



evaporation tim escale is given by (|Lubow & Shu||1975 
Lai et al.||2010D 

GM^ /ON 

'TRoche ~ 5-^, (2j 

where Prhui is the atmospheric density at the Hill ra- 
dius and c is the speed of sound. Generally, mass-loss by 
Roche-lobe overfiow is orders of magnitude faster than 
that by extreme ultraviolet radiation-driven winds. Al- 
though any criterion based on the above estimates is only 
accurate to within an order of magnitude or so, we find 
this to be sufficient for our purposes, as we typically find 
a rapid transition from r ^ Gyr to r <C Gyr across two 
models that neighbor each-other in mass. 

The mass-radius relationships for planets with core 
masses of 1,3 and SMq are presented in panels A, B, 
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and C of Figure ([3| respectively. Black dots represent the 
radii obtained through numerical experiments while the 
curves depict interpolation functions that run through 
the data. Thick curves imply models that are secure 
against evaporation while thin lines depict unstable mod- 
els. In addition to the irradiated models (shown with 
blue lines), isolated (i.e. no irradiation) models are also 
presented and are shown with black lines. For reference, 
Jupiter's, Saturn's and Neptune's radii are also marked. 

The results highlight the fact that accounting for stel- 
lar irradiation, giant planetary radii can persist to sur- 
prisingly low masses (that is, M < IOMq). Figure ([3| 
further affirms that the behavior of planetary structure 
is largely dictated by the associated core mass. Note 
that all models with a SMq core are stable against evap- 
oration and roughly follow the cold (i.e. isolated) mass- 
radius relationship. On the contrary, IM^ core mod- 
els are largely unstable below M < ISMq but can have 
radii comparable to that of Jupiter prior to the onset 
of evaporation. A similar scenario is observed for the 
3M0, Teq = 500K set of models. Indeed, these models 
are essentially always characterized hy R '^ ^jup above 
M>8Me. 

It is interesting to note that some of our models (e.g. 
those corresponding to 1 and SMq and Teq = TOOK) have 
radii that are bigger than that of Jupiter. As already dis- 
cussed above, this upturn in radii is a direct consequence 
of the softening of the equation of state at lower pres- 
sures (an ideal gas has a softer equation of state than 
the deep interior of Jupit er). While reminiscent of the 
inflated Hot J upiter radii ( |Guillot|[2QQ5l [Fortney fc Net- 
telmann|2QlQ ), these objects are tundanientally ditterent, 
since they require no additional heat sources or mecha- 
nisms for stalling gravitational contraction. That said, 
it is unclear if such objects are particularly significant 
within the context of the observational sample, since the 
models that show such an excess are close to the evapo- 
ration boundaries of the mass-radius diagrams. In fact, 
accounting for coupled evolution of gravitational contrac- 
tion and mass-loss in a more self-consistent matter will 
likely yield an exclusion region that is a bit larger than 
what is shown in Figure 3. 

4. DISCUSSION 

In this letter, we have examined the structure of mod- 
erately irradiated low-mass low-density extrasolar plan- 
ets. We began by analyzing the interior of a compara- 
tively well characterized planet iTep/er-SOd, and showed 
that the planet is likely composed of an extensive gaseous 
H/He envelope, surrounding a core that makes up less 
than a third of its total mass. Although qualitatively 
this object resembles a scaled down version of Saturn, 
it is important to recall that the mass of i^ep/er-30d is 
typical of much more metal-rich objects such as Uranus 
or Neptune. The existence of i^ep/er-30d immediately 
suggests that range of planetary interior configurations 
that occur in nature is much wider than that available 
for detailed study within the realm of the solar system. 

Prompted by this notion, we extended our calculations 
to quantify planetary mass-radius relationships for cored 
low-mass gaseous planetary objects of solar composition. 
Our calculations underline the importance of stellar ir- 
radiation on the evolutionary tracks of low-mass objects. 
In particular, the constructed mass-radius relationships 
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Fig. 3. — Mass- Radius relationships of low-mass, gas-dominated 
planetary models. Panels A, B, and C correspond to planets with 
core-masses of Mcore = 1,3 and SMq respectively. On each panel, 
mass-radius relationships corresponding to equilibrium irradiation 
temperatures of Teq = 300, 500 and TOOK are shown as blue lines. 
Additionally, isolated mass-radius relationships are shown as black 
lines. Solid lines run through models that are stable against evap- 
oration while the converse is true for thin lines. Note that radii 
characteristic of giant planets are readily attainable for mildly ir- 
radiated M ~ IOM0, Mcore = 1, 3M0 planets. 



suggest that the radius of an irradiated body may exceed 
that of its isolated counterpart by as much as a factor 
of -- 2 (e.g. the case of M :^ lOM®, Mcore = SM®, 
Teq = 500), bringing the radius well into the character- 
istic giant planet range. Collectively, our results suggest 
that extreme care must be taken in the interpretation 
of giant transit radii from the Kepler sample, since the 
mass range corresponding to such radii can be quite ex- 
tensive (i.e. spanning almost two orders of magnitude). 
One may wish to argue against a significant popula- 
tion of bodies like those considered in this work based on 
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the (im) probability of their formation, since the gaseous 
component of our models is much enhanced over the stan- 
dard models of typical objects in the considered mass 
range. Indeed, it is often said that one must have a 
"critical" core mass of M^ove ^ IOMq in order to trigger 
gas accretion. However, this claim is ill-founded and is 
not actually relevant since a hydrostatically supported 
atmosphere around a core can be more massive than 
that envisioned within the context of the standard mod- 
els (Poll ack et al.||199"6 ) if either accretion is slower, the 
molecular weight ot the envelope is larger, or the opacity 
is increased. This is evident for example in the simple 



analytical mode l s of Ste venson (1982b) (see also Ikoma 
fc Genda| (|2QQ6[); |Broeg| (2009)). Protoplanets may also 



have circumplanetary disks that qualitatively change the 
characteristic accretion pattern and affect the planetary 
energy loss. Furthermore, alternative formation scenar- 
ios could likely be envisioned, a speculative example be- 
ing one where objects of this type are sculpted out of 
more massive planets by intense ultraviolet-driven mass- 
loss during the first ~ lOOMyr of the stellar lifetime. In- 
deed, such scenarios have alr eady been proposed in the 



exo-p lanetary context (Baraffe et al.||2006[ [Lopez et al. 
2012J ) 



Ultimately, our aim here is not to argue for or against 
any particular formation scenario for sub-Neptune 
mass gas-dominated planets. Rather, similarly to what 
has been done for iTep/er-SOd, we propose that their 
existence can be validated or ruled out observationally. 
Beyond standard methods like transit timing variations. 



the most obvious approach to this is through radial- 
velocity monitoring of transiting planets. That is, if a 
giant-plenet-like radius is firmly established for a given 
object through transit observations but a commensurate 
radial velocity signal is not observed in the host star, 
such an object is likely characterized by a very low mass. 
Another approach to mass discrimination is exclusively 
photometric, and takes advantage of dependence of 
the transit radius on spectral frequency. Although 
the atmospheric scale-heights of Hot Jupiters comprise 
~ 1% of their radii at most, for ~ IOMq (albeit a factor 
of ~ 3 cooler) planets with similar radii, the scale height 
is increased by about an order of magnitude. As a 
result, the chord optical depth of unity may correspond 
to substantially different radii in visible and infrared 
light. Both of these observational avenues should 
become readily available as the radial velocity precision 
continues to improve and future space-based missions 
such as JWST commence. 
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